Spin-Hall effect and spin-Coulomb drag in doped semiconductors 
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In this review, we describe in detail two important spin-transport phenomena: the extrinsic 
spin-Hall effect (coming from spin-orbit interactions between electrons and impurities) and the 
spin-Coulomb drag. The interplay of these two phenomena is analyzed. In particular, we discuss 
the influence of scattering between electrons with opposite spins on the spin current and the spin 
accumulation produced by the spin-Hall effect. Future challenges and open questions are briefly 
discussed. 



I. INTRODUCTION. 

In recent years two important spin-transport 
phenomenal have been discovered in semiconduc- 
tors, in the conducting (metallic) regime the spin-Hall 
effect (SHE)and the spin-Coulomb drag. 

The spin-Hall effect^ is a bulk property of the 
semiconductor with a strong spin-orbit interactions 
in the metallic regime. SHE is a close cousin of 
the anomalous Hall effect. Anomalous Hall effect 
(AHEj^MI&a^li^lllilSiia . g ^ generation of a 

transverse charge and spin polarization current in re- 
sponse to an electric field. It appears in ferromagnets 
with strong spin-orbit interactions like GaMnAs. In con- 
trast, the spin-Hall effect (SHE) is the generation of a 
transverse spin polarization current alone in response to 
an electric field in a paramagnetic medium with spin- 
orbit interactions and in the absence of a magnetic field. 
By analogy with AHE there are two mechanisms gen- 
erating SHE: impurities and band structure. While the 
impurity mechanism was suggested many years ago by 
Dyakonov and Pere l 17 ' 18 ' 19 ' 20 , the second mechanism, 
originating from band structure, has come to the fore- 
front in recent year o 21 i 22 . Again, similar to AHE, a lively 
debate arose about which of these two mechanisms - 
extrinsic (coming from impurities) or intrinsic (coming 
from band structure) 2 ^^^^* 2 ^^^ 3 -^^ - 
is more important, and how to distinguish experimentally 
between the two 37 i 38 . Because of the transverse charge 
response that comes with it, the AHE can be detected by 
purely electrical measurements. However, this is not the 
case for the SHE, because the spin-polarization current 
can not be directly measured in transport. The spin-Hall 
effect in semiconductors (GaAs, ZnSe) has been mainly 
observed in optical experiment s 39 ' 40 ' 41 ' 42 ' 43 . The idea of 
these experiments is following: in a finite size sample, 
charge current induces a transverse uniform gradient of 
spin density (spin accumulation) which increases until 
the steady-state is achieved. This spin accumulation can 
be measured quite clearly by observing a change in po- 
larization of a reflected beam of light (Kerr effect). This 
method has been successfully applied to n-type GaAs 



sample s 39 ' i ' . In another experiment performed on 
p-type GaAs, the spin-accumulation was revealed by the 
polarization of the recombination radiation of electrons 
and holes in a two-dimensional LED structure. More re- 
cently, it has become possible to study the time evolution 
of optically injected charge and spin currents 4 ^, and to 
monitor the dynamics of spin accumulation in semicon- 
ductors^ 3 -. 

The possibility of detecting the SHE by electrical mea- 
surements in mesoscopic samples was theoretically sug- 
gested i n 19 i 45 . In that proposal, an electric current driven 
in one of the legs of an H-shaped structure generates a 
transverse spin current in the connecting part due to the 
SHE. Then, due to inverse spin-Hall effect, this spin cur- 
rent produces a voltage difference in the second leg of 
the structure 4 ^. Very recently this proposal was real- 
ized in the H-shaped structures of the size of one mi- 
crometer, fabricated on the HgTe/CdTe quantum wells 
in the inverted regime 4 ^. These quantum wells charac- 
terize a very long mean free path (larger than a couple 
of micrometer) and a very strong spin-orbit coupling 47 
of the intrinsic (Rashba^ 8 - ) type and authors concluded 
that observed voltage (of the order of microvolt) is the 
proof of the first measurement of ballistic spin-Hall ef- 
fect in transport. In a similar setup but using the in- 
verse spin-Hall effect alone* 9 -, the spin-polarization was 
converted to electric signal and at least an order of mag- 
nitude weaker electrical signal was detected in metals, 
such as A l 50 i 51 . The origin of spin-Hall effect in metals 
is still under debat o 52 ' 53 ' 54 . Although theoretically esti- 
mated intrinsic contribution in some of these materials 
(like platinum) can be larg o 53 ' 54 , the calculations were 
performed, so far, for macroscopic systems and did not 
include an extrinsic contribution. Therefore further ef- 
fort is warranted in this field, aimed at clear-cut distinc- 
tion between different mechanisms contributing to the 
total spin-Hall effect in metals and semiconductor s 37 ' 38 . 
In this review we mainly focus on the extrinsic spin-Hall 
effect. 

The quantum spin Hall (QSH) state is a novel topolog- 
ically nontrivial insulating state in semiconductors with 
strong spin-orbit interaction o 55 ' 56 ' 57 ' 58 ' 59 ' 60 ' 61 , very dif- 
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ferent from the SHE. The QSH state, similar to the 
quantum Hall (QH) state, has a charge excitation gap 
in the bulk. However, in contrast to the QH state, the 
QSH state does not require existence of the magnetic 
field. Therefore for the QSH state, time reversal sym- 
metry is not broken and instead of one spin degenerate 
edge channel (as in the QH effect), two states with oppo- 
site spin-polarization counterpropagate at a given edge. 
The QSH effect was first proposed by Kane and Mele 
for graphene^. However, the gap opened by the spin- 
orbit interaction turned out to be extremely small on the 
order of 10~ 3 meV. Very recently Bernevig and Zhang 
theoretically proposed that the QSH effect should be vis- 
ible in HgTc/CdTe quantum wells with inverted band 
structure^ 7 - and the experimental discovery of QSH ef- 
fect in this material followed shortly afterwards^ 8 . QSH 
effect and the SHE effect are two distinct phenomena. 
While transport in the QSH effect occurs in the spin 
edge channels of an insulating material, the SHE involves 
transport in the bulk of a conductor. This review is fo- 
cused on the semiconductor spin transport in the metal- 
lic regime, where the bulk is conducting. Specifically, we 
will summarize here the current status of our knowledge 
concerning two important spin transport phenomena in 
this regime: the spin-Hall effect and the spin-Coulomb 
drag. We refer a reader to the recent review^ for further 
details concerning the QSH effect. 

The spin-Coulomb dra$2£2£2£Z£ML£2. is a many- 
body effect arising from the interaction between electrons 
with opposite spins, which tends to suppress the relative 
motion of electrons with different spins and thus to re- 
duce the spin diffusion constant. This effect has been 
recently observed in a (110) GaAs quantum well (which 
is essentially free of spin-orbit interaction) by Weber et 
alJ£ by monitoring the time evolution of a spatially vary- 
ing pattern of spin polarization, i.e. a spin grating. The 
rate of decay of the amplitude grows in proportion to the 
square of the wave vector of the grating, and the coeffi- 
cient of proportionality is just the spin diffusion constant. 
The measured value of the spin diffusion constant turns 
out to be much smaller than the single particle diffu- 
sion constant (deduced from the electrical mobility) and 
the difference can be quantitatively explained in terms 
of Coulomb scattering between electrons of opposite spin 
orientation drifting in opposite direction, thus lending 
support to the theory of spin-Coulomb drag, as described 
in detail in Section III. 

Intuitively when both the spin-Hall effect and the spin- 
Coulomb drag are present, the spin current generated by 
SHE should be reduced and therefore it is important to 
take a look at the combined influence of these two effects 
on spin transport. 

The rest of the review is organized as follows: in Sec- 
tion II we describe the extrinsic spin-Hall effect using the 
Boltzmann equation approach; in Section III we describe 
the spin-Coulomb drag effect; in Section IV we discuss 
in detail the influence of spin-Coulomb drag on the ex- 
trinsic spin-Hall effect; in Section V we briefly describe 



the intrinsic spin-Hall effect and the influence of intrinsic 
spin-orbit coupling on the spin-Coulomb drag; in Sec- 
tion VI, we discuss possible scenarios for the evolution 
of the SHE in semiconductors as a function of mobility. 
Conclusions and open challenges are presented in Section 
VII. 



II. EXTRINSIC SPIN-HALL EFFECT. 

Fig. (1) shows a setup for the measurement of the SHE. 
An electric field (in the x direction) is applied to a non- 
magnetic two dimensional electron gas (2DEG). In re- 
sponse to this, a spin current begins to flow in a direction 
perpendicular to the electric field (the y direction, see 
Fig. (1)): that is to say, spin up and spin down electrons 
with "up" and "down" defined in respect to the normal 
to the plane, drift in opposite directions perpendicular to 
the electric field. 
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FIG. 1: Schematics of the spin-Hall effect. An electric field 
applied along the x axis, E x , induces a transverse spin current 
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There are two extrinsic mechanisms of generation 
of transverse spin current: skew scattering and side 
jump^- i 70 ' 71 . They both arise from the effect of the spin- 
orbit interaction on electron-impurity collisions. Skew 
-scattering arises from the asymmetry of the electron- 
impurity scattering in the presence of spin-orbit interac- 
tions: electrons that are drifting in +x direction under 
the action of electric field are more likely to be scattered 
to the left than to the right if, say their spin is up, while 
the reverse is true if their spin is down. This generates 
a net z-spin current in the y direction. This mechanism 
is also known as "Mott scattering"— and has been long 
known as a method to produce spin-polarized beams of 
particles. 

The second effect is more subtle and is caused by the 
anomalous relationship between the physical and canon- 
ical position operator, as will be explained below. It is 
called "side-jump" , because semi-classically it can be de- 
rived from a lateral shift in the position of wave packet 
during collision with impurities. Without resorting to 
this description, we could arrive at the same side-jump 
term starting from the quantum kinetic equation and in- 
cluding the anomalous part of position operator. Figs. 2 
and 3 present simple pictures of skew scattering and side 
jump mechanisms respectively. 
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operator is modified, taking the form 
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FIG. 2: Schematic picture of the skew-scattering effect. Due 
to the spin-orbit interaction between electrons and impuri- 
ties, electrons with different spin orientations are deflected to 
opposite edges of sample. 
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FIG. 3: Schematic picture of the side-jump in a head-on colli- 
sion with the impurity. The center of electron wave packet is 
shifted by Ar p h ys in a direction perpendicular to the change 
of momentum during the collision. The physical origin of this 
shift is described in the text. 



After this brief, pictorial presentation of different 
mechanisms contributing to the extrinsic spin-Hall effect, 
we are now ready to begin a more detailed analysis. 

First of all, because both mechanisms depend cru- 
cially on the spin-orbit interaction, it is necessary to 
say something about the character of this interaction in 
the solid state. The derivation of this interaction in- 
volves steps that closely parallel the derivation of the 
spin-orbit interaction from the Dirac equation for a sin- 
gle electron in vacuum. In that case, we arrive at the 
effective one band (Pauli) Hamiltonian by applying a uni- 
tary transformation^ that decouples the electrons from 
the positrons, and then projecting onto the electron sub- 
space. However, during this transformation the position 



where aoh = h 2 /Am 2 c 2 w —3.7 x 10~ 6 A 2 is the strength 
of spin-orbit interaction for bare electron in vacuum and 
m e is the bare electron mass. This physical position op- 
erator is valid for electrons in a conduction band. For a 
general case, one has to apply the form of r p h ys as shown 
in Eq. (74) of paper—. Substitution of the modified posi- 
tion operator in the potential, followed by an expansion 
to first order in ao, leads to the standard form of the 
spin-orbit interaction in vacuum: 



Hso = a (px W(f)) 



(2) 



where V(r) is the external potential acting on an elec- 
tron, p is the electron momentum and a is the vector 
of the Pauli matrices. In semiconductors, the spin-orbit 
(SO) interactions play double role. First, we have SO 
effects induced by the periodic crystal potential Vo(r). 
This causes splitting of the p-like valence band at k = 0, 
in semiconductors like GaAs, into a fourfold degenerate 
band with total angular momentum j = 3/2 (heavy and 
light hole bands) and "split-off band" with j = 1/2. Fur- 
ther, the periodic crystal field gives rise to a small spin- 



orbit interaction of the order of 



4m^c 2 



[px W (M)] 



on electrons in the conduction band. Second, there is the 
SO interaction induced by any external potential (differ- 
ent than Vo(r))if one wants to find an effective model say 
for the conduction band. In other words, if we perform 
a similar unitary transformation (as we did for the Dirac 
model) in a semiconductor, folding eight bands (conduc- 
tion band, heavy and light hole bands as well as split-off 
valence band) into an effective model for the conduc- 
tion band then the resulting spin-orbit interaction will 
be again of the same form as in Eq. (T2J) , but with a much 
larger "coupling constant" >2£ 
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hP 2 

3m 2 
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El {E g + Aso) 2 



(3) 



where E g is the gap energy between conduction and 
heavy/light holes bands, A50 is the splitting energy be- 
tween heavy/light holes and split-off bands, P is the ma- 
trix element of the momentum operator between the con- 
duction and the valence-band edges. Using values of the 
parameters appropriate for the 2DEG in Alo.iGao.gAs^S 
with a conduction band mass m — 0.074m e we find 
ah = AAA 2 . Therefore in the conduction band of semi- 
conductors the spin-orbit interaction is six orders of mag- 
nitude larger than in vacuum and has the opposite sign. 
Obviously the spin-orbit interaction induced by the peri- 
odic crystal field on the conduction band can be omitted 
as a correction many orders of magnitude smaller. How- 
ever, ah is much smaller than the square of the effec- 
tive Bohr radius in GaAs ( 10 4 A 2 ) and in this sense the 
spin-orbit coupling can still be considered a small per- 
turbation. Notice that the form of the physical position 
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operator for the conduction band is described by Eq. JTJ) , 
with a replaced by a. 

Taking this into account and omitting for the time be- 
ing the electron-electron interaction, we see that our ef- 
fective Hamiltonian for electrons in the conduction band 
of GaAs takes the form 



H ni — H + V e i(r p hys) + E(r p h V s), 



where 



(4) 



(5) 



is the kinetic energy of electrons in conduction band (m 
being the effective mass of the conduction band) and 



Vei(r p hys) ^ V ei (n) + a ^ \pj X ViV e i(ri) 



(6) 



where V e i(?i) is the impurity potential and the spin-orbit 
interaction between electrons and impurities arises from 
the Taylor expansion of the impurity potential around 
the canonical position operator r*j. Finally, 



where we took the interaction with the external electric 
field E to be eE ■ r p h ys ■ 

The physical velocity operator is the time derivative of 
the physical position operator i.e.: 



Vphys 

and has the form 



Vphys i HI 



(8) 
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ViVeiifi) + eE 



dpi , n v 
x o-j - a— x a h (9) 
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where the first two terms on the right hand side are 
derivatives of the canonical position operator while the 
last term originates from the time derivative of the 
anomalous part of position operator. However, since 
the total force Fi — dpi/dt consists of a force origi- 
nating from impurities and one from the electric field 
Fi = — WiV e i(fi) — eE, the second and last terms of 
Eq. ([9]) are equivalent. Therefore v p h ys can be written 
in the following compact form: 
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2a 



ViV ei (n) + eE 



X CTj 



(10) 



One can see that in our model the z-component of spin 
is conserved because it commutes with the Hamiltonian. 
We exploit the conservation of a z by defining the quasi- 
classical one-particle distribution function f a (r, k,t), i.e. 
the probability of finding an electron with z-component 
of the spin S z — ^cr, with a = ±1, at position r with 

momentum p = hk at the time t. In this review we 



focus on spatially homogeneous steady-state situations, 
in which f a does not depend on r and t. We write 

U(f,k,t) = foa(e k ) + fiAk) , (11) 
where /oo-(efc) is the equilibrium distribution function - a 
function of the free particle energy e k = - and f\ a (k) 
is a small deviation from equilibrium induced by the ap- 
plication of steady electric fields E G (a = ±1) which cou- 
ple independently to each of the two spin components. In 
the next few sections we will apply the Boltzmann equa- 
tion approach to calculate f± taking spin Hall effect and 
spin Coulomb drag into account on equal footing^ 

To first order in E a the Boltzmann equation takes the 
form 



hk 



eE a fd a (e k ) = fia(k) c 



(12) 



where /^(efc) is the first derivative of equilibrium distri- 
bution function with respect to the energy and fia(k) c is 
the first-order in E a part of the collisional time derivative 
f a {k) c due to various scattering mechanisms. For the 
electron-impurity scattering mechanism the collisional 
time derivative has the following form: 



E{r vhya )~^leE-?i+ea($iXE)-&i\ (7) f a [k) c 



where W^p is the scattering rate for a spin-er electron 

to go from k to k' , and e^a is the particle energy, includ- 
ing an additional spin-orbit interaction energy due the 
electric field: 



ffcer = (k + 2eaha(E a X z) ■ k , 



(14) 



The peculiar form of t ka , which differs from the naive 
expectation e k + eE ■ f p h ys by a factor 2 in the second 
term, is absolutely vital for a correct treatment of the 
"side-jump" contribution. The reason for the factor 2 is 
that the 6 function in Eq. (|13p expresses the conservation 
of energy in a scattering process. Scattering is a time- 
dependent process: therefore the correct expression for 
the change in position of the electron Ar p h ys must be 
calculated as the integral of the velocity over time: 



Arphys 



V p liysdt 



(15) 



Before solving integral in Eq. (|15p. let's think for a mo- 
ment about scattering event. Let's take Ap = p out — p- m 
(see Fig. (3)) to be the change in momentum of an elec- 
tron wave packet during collision with an impurity. Dur- 
ing the very short time of collision, W(r) = —dp/dt is 
very large and therefore second term of Eq. pTJ|) com- 
pletely dominates the velocity. Therefore, we disregard 
first term in the velocity formula and obtain a following 
form for the electron wave packet displacement: 



Af ph y S = -2aAp x a/fi 
Therefore, Eqs. (|9I14I15[) are consistent. 
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A. Skew-scattering 



From the general scattering theory, developed for in- 
stance ir>Z£, one can deduce the form of scattering prob- 



ability from k to k 



m,76 



a.s: 



and 



1 T->o mA f 2n 

IT ~ ^ Jn 



, d9 W a {k F ,6) sin 2 9 (22) 



In Fig. (4) and Fig. (5), we present the integrated sym- 



W tk> + uW U> (* x *')* - ) > (17) 



kk 



where for centrally symmetric scattering potentials Wjip 

and WSp depend on the magnitude of vectors k and 

kl and the angle 9 between them. Furthermore the 
left/right asymmetry of skew scattering is included ex- 
plicitly in the factor (k x k') z and therefore both Wi~ 

and W2p are symmetric under interchange of k and k'. 
Taking into account the form of the scattering probability 
as well as the conservation energy during the scattering 
process we obtain following expression for the linearized 
collisional derivative: 



c,zmp 



--J2 w kk 
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Wlp{fUk)-fiAk')}5(e k -e k ,) 



k,(k x k% {/ lCT (fe) + f la (k')}s(e k - ey) 



-2cr 



k' 



Wl p f^(e k )ean(E a x z) ■ (k - k')S(e k - e k> ) 



(18) 



where fi a (k) = — f' 0(y (t k )Tik ■ V a (k) and the first term on 
r.h.s. of this formula is the symmetric scattering term, 
the second one is the skew scattering term, while the last 
one will be ultimately responsible for the side jump. To 
find the drift velocity, V a {k), we need to multiply both 
sides of Boltzmann equation (|12p by hk/m and integrate 
over k space, and therefore derive V a (k) self-consistently 
from the condition: 



Ehk 
— 
m 



m 



foa-( 6 k) — ^ t fla(k) c ,imp ■ 



(19) 

After substituting of collisional derivative Eq. (TT8|) into 
Eq. (|19p and with the assumption that we can omit the 
k-dependence of drift velocity in low temperatures, we 
arrive at the following formula for V a to first order in the 
spin-orbit interaction: 



E a — a—E a x z 



2eaa{E a x z) , (20) 



where in the limit of zero temperature the symmetric 
scattering rate 1/r and the skew scattering rate 1/t ss 
simplify to: 



1 t-^o mA 



4ir 2 h 2 



d9 W s (k F ,9)(l- cos 9) , (21) 
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FIG. 4: Integrated symmetric scattering rate in units of 
2mh/m 2 A as a function of k 2 for a model circular well at- 
tractive potential Vo = — 5meV and radius a — 9.45nm (de- 
scribed in Reference^). We choose the parameters typical 
for the experimental 2DEG confined in Alo.1Gao.9As quan- 
tum well i.e. density of electrons and impurities U2d = rii = 
2.0 x 10 12 cm" 2 , m = 0.074m e , and mobility fx = 0.1m 2 /V.s. 

o 2 

The effective spin-orbit coupling ah = 4AA in accordance 
withZi. 

metric and asymmetric scattering rates calculated from 
Eqs. pTjl and (|22[) for a simple step potential of the form: 



V(r) = V Q 9{a - r) + aaL z S z S(r - a)V 



(23) 



where Vq is the attractive electron-impurity potential, 
a = ah I a 2 is the renormalized spin-orbit interaction, a 
is the effective spin-orbit coupling constant for the con- 
duction band, a is the impurity radius, L z and S z are 
the orbital angular and spin angular momenta, respec- 
tively. The simple step potential described in Eq. (|23|) 
is presented in Fig. (6). For parameters typical of III-V 
semiconductors, the asymmetric and symmetric scatter- 
ing rates are almost flat as functions of the energy of the 
incoming electron (see Figs. (4,5)). Comparing symmet- 
ric and asymmetric scattering rates one can find directly 
from picture, that for typical densities of 2DEG and typ- 
ical Bohr radii in semiconductors the ratio of symmetric 
scattering time to skew scattering time is: 



0.002-0.003 



(24) 



We assumed, in above discussion, that the electron- 
impurity potential is attractive, which is the most com- 
mon case in semiconductors, where ionized donors play 
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FIG. 5: Integrated asymmetric scattering rate in units of 
2mh/m 2 A as a function of k 2 for a model circular well at- 
tractive potential Vo = — 5meV and radius a = 9.45nm (de- 
scribed in Reference^). We choose the parameters typical 
for the experimental 2DEG confined in Alo.1Gao.9As quan- 
tum well i.e. density of electrons and impurities 7120 = rii = 
2.0 x 10 12 cm" 2 , m = 0.074m e , and mobility fj, = 0.1m 2 /V.s. 
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FIG. 6: The electron-impurity potential V(r) used to estimate 
the ratio of symmetric to asymmetric scattering. 



the role of impurities. However, the sign and the ampli- 
tude of skew scattering contribution (and more precisely 
skew scattering rate) depends strongly on the electron- 
impurity potential. 



B. Side jump current and resistivity matrix 

Now to determine the side jump contribution to the 
spin-Hall effect we need to carefully define the spin- 
current. The spin-current density operator must be cal- 
culated taking into account the form of the physical ve- 
locity (see Eq. (|TTJ|l V 
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where the factor 2aFi X /h vanishes because the net force 
Ft acting on an electron is zero when averaged over a 
steady-state ensemble. Therefore, the spin component of 
the current is 



Jo 
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(26) 



and using Eq. (|20|) for the drift velocity we obtain the 
following form of side-jump current: 
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2e 2 aan a E a 



(27) 



which evidently arises from the last term of Eq. (|20[) . 
This simple form of the side-jump current is valid only for 
electrons in the conduction band. In this case, the side- 
jump current depends only on the spin-orbit coupling, 
the density of electrons and the spin-dependent electric 
field. The complete relation between the spin-component 
of current (from side jump and skew scattering contribu- 
tions) and electric field has the following form: 



x z 



(28) 



where p l a 



is the Drude resistivity, p s a 



is the skew scattering resistivity, and the last term in 
square brackets is the side-jump contribution to the re- 
sistivity: X a = 22121. Eq. $ZE§ yields to the following 
resistivity tensor (in the basis X|, j/j, X[, y^): 



P 





Vo 



„SS 1 \ „i 

-p t + A T p T 



D 



p 

p\ 






Pf 



-pT + hPi 
Pt s -hp? p? 

(29) 

The diagonal part of the resistivity reduces to the Drude 
formula p^ — ™ , as expected. The spin-orbit in- 
teraction is entirely responsible for the appearance of an 
off-diagonal (transverse) resistivity. The latter consists of 
two competing terms associated with side-jump (A a p®) 
and skew-scattering (p* s ), as seen in Eq. (p9|) . The signs 
of side-jump and skew-scattering terms are opposite for 
attractive electron-impurity potential. Although, this is 
the most typical case in doped semiconductors, it is im- 
portant to emphasize that side-jump and skew scattering 
terms have equal signs in the case of a repulsive electron- 
impurity potential. Also, as one can see from tensor |29|) 
contributions scale differently with the mobility. Since 
t ss ~ r the skew scattering contribution to the resistiv- 
ity is proportional to where p is a mobility, while the 
side jump contribution scales as 1/p 2 - The opposite signs 
of two contributions, and different scaling with mobility 
could allow to distinguish between them in the experi- 
ments (see Section IV)2£. As expected, in the absence of 
electron-electron interactions, the resistivity tensor 
does not include elements between the opposite spins. 
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III. SPIN-COULOMB DRAG 

Ordinary Coulomb drag is caused by momentum ex- 
change between electrons residing in two separate 2D lay- 
ers and interacting via the Coulomb interaction (for re- 
view see2£). The Spin-Coulomb drag is the single-layer 
analogue of the ordinary Coulomb drag. In this case spin- 
up and spin-down electrons play the role of electrons in 
different layers and the friction arises (due to Coulomb 
interactions) when spin-up and spin-down electrons move 
within one single layer with different drift velocities^. 

The simplest description of the spin-Coulomb drag is 
given in terms of a phenomenological friction coefficient 
7. Later in the section we will show that the Boltz- 
mann equation approach confirms this phenomenological 
description. Let us start with the equation of motion for 
the drift velocity of spin-cr electrons: 



where p SD = mj/ne 2 is the spin Coulomb drag resistiv- 
ity and p' a = m/n a e 2 T' a . Several features of this matrix 
are noteworthy. First the matrix is symmetric. Second 
the off-diagonal terms are negative. The minus sign can 
be easily explained, p^ is the electric field induced in 
the up-spin channel by a current flowing in the down- 
spin channel when the up spin current is zero. Since a 
down spin current in the positive direction tends to drag 
along the up-spins, a negative electric field is needed to 
maintain the zero value of the up-spin current. There is 
no limit on the magnitude of psd- The only restriction 
is that the eigenvalues of the real part of the resistivity 
matrix should be positive to ensure positivity of dissi- 
pation. Finally, the spin-Coulomb drag appears in both 
diagonal and off-diagonal terms so the total contribution 
cancels to zero (in accordance with Eq. (j3"Tj) ') if the drift 
velocities of up and down spins are equal. 



mV 



-eN a E a + F a 



(30) 



where N a is the number of electrons with spin a, F a ^ a 
is the net force exerted by — a spins on a spins, — is the 
rate of change of momentum of electrons with spin a due 
to electron-impurity scattering and is basically the Drude 

scattering rate, -^f- is the rate of change of momentum 
due to electron-impurity scattering in which electron flip 
its spin from — a to a. From Newton's third law one 
immediately sees: 



(31) 



and by Galilean invariance this force can only depend on 
the relative velocity of the two components. Hence, for 
a weak Coulomb coupling one writes: 



F a -a = -lN a (Vi, a - Vi-a) - 



(32) 



Coulomb scattering 



Spin-flip scattering 



p+qt 





1/T' 



FIG. 7: Comparison of Coulomb scattering with the spin-flip 
scattering. At finite temperature Coulomb scattering can be 
a more effective mechanism of momentum exchange between 
up- and down-spin populations than spin-flip collisions with 
impurities. 



where n a is the density of electrons with spin a and 7 is 
a spin-drag scattering rate. Taking into account Eq. (j2l)|) 
and applying Fourier transformation to Eq. (|30[) one gets 
following equation on the spin current density: 

iujj a {uj) = + 7 + — j CT (w) 

m \ n T a J 

7— + V)7-"M ( 33 ) 

n t' 



Inverting Eq (|33|) gives us electric field: 



E a [u) 



m rrv) ri- a 



n a e 2 n a e 2 T a ne 2 n 



m m7 
7 + — 2 



n a e 2 T' (J in- 



(34) 



From this we can immediately read the resistivity tensor. 
Its real part, in the basis of x-\, X[, has following form: 

(pf+P^n./m -p sd - P \ \ 
P ~\ -p SD -p\ pf+p^/nj (35) 



Let us take a closer look at the competing off-diagonal 
terms: spin-Coulomb drag and spin-flip resistivities. At 
very low temperature spin flip processes win because 
in this limit the Coulomb scattering is suppressed by 
phase space restrictions (Pauli's exclusion principle) and 
7 tends to zero as T 2 in three dimensions and T 2 lnT 
in 2D. However, the spin-flip processes from electron- 
impurity collisions do not effectively contribute to mo- 
mentum transfer between two spin-channels. An up-spin 
electron that collides with an impurity and flips its spin 
orientation from up to down is almost equally likely to 
emerge in any direction, as shown in Fig. (7), so the mo- 
mentum transfer from the up to the down spin orienta- 
tion is minimal and independent of what the down spins 
are doing. However, the situation looks quite different 
for electron-electron collisions: the collision of an up-spin 
electron with a down-spin electron leads to a momentum 
transfer that is preferentially oriented against the relative 
velocity of the two electrons (see Fig. (7)) and is propor- 
tional to the latter. Taking the spin-flip relaxation time 
to be of the order of 500ns (ten times larger than spin- 
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relaxation time in GaA a 79 i 80 ), and the value of - of the 

order of the Drude scattering time, around 1 p o 65 ' 66 , and 
temperatures of the order of the Fermi energy Tp ~ 3001^ 
we estimate that the spin-Coulomb drag contribution will 
dominate already for T ~ 0.3K. Further, for mobilities 
typical for semiconductors 10 4 ~ 10 5 cm 2 /Vs the ratio 



of spin-Coulomb drag to the Drude resistivity can be as 
large as ten (see more detailed discussion in the next Sec- 
tion and Fig. (9)). 

Returning to Boltzmann approach, the electron- 
electron contribution to the collisional derivative has the 
form 65 



U(k)c,e- e --J2 W c (ka,p- a;k'aj - a) [1 - f a (k')][l - f_„tf)\ 

k'pp' 

-f*(k')f-«(?)[l - /.(fc)][l - f-AP)}} k+pfi+pS&tor + ep-. ~ *** ~ tp'-a) , (36) 

I 



where Wc(ka,p — a; k'a^ff — a) is the electron-electron 
scattering rate from fee, p — a to k'a, p 1 — a, and the Pauli 
factors f a (k), 1 — f a (k') etc. ensure that the initial states 
are occupied and the final states empty as required by 
Pauli's exclusion principle. Notice that, for our purposes, 
only collisions between electrons of opposite spins are rel- 



evant, since collision between same-spin electrons con- 
serve the total momentum of each spin component. Af- 
ter substituting the linearized Boltzmann equation into 
Eq. (|36p and in the absence of spin-orbit interactions one 
derives the following Coulomb collision integral 4 ^: 



U(k) c<e - e -—T^Yl W c (ka,p- a; k'a ;0 - a)[hV a - W-^k-k 1 ) 



k B T , 

k'pp' 

/o(r(efe)/o-cr(ep)/o(T(— efc')/o-cr(— e p')H+p,k'+p'^ eka + e P~°" _ ek ' a ~ e p'~^ 

I 



(37) 



where T is the temperature, k B is the Boltzmann 
constant, and we have made use of the identity 

/(kr(efc)/0-a(ep)[l - /0*M][1 - /0-«r(v)] = I 1 : 
/0<r(£fc)][l - fo-*(£p)}foa(tk')fo-*(£p') for ^ka + £p-a - 



tk'a — e p '-a — 0. The collision integral Eq (|37j) is pro- 
portional to the difference of velocities for spin-up and 
spin-down electrons. Therefore if a finite spin current is 
set up through the application of an external field, then 
the Coulomb interaction will tend to equalize the net 
momenta of the two spin components, causing (Vj-) -(Vj,) 
to decay and thus can be interpreted, as we explained 
before, as a damping mechanism for spin-current. 

After substituting Eq. (f3"T)) into self-consistent equa- 
tion for a drift velocity Eq. (fT9|) . one obtains the spin- 
drag coefficient 7 i.e. the rate of momentum transfer 
between up and down spin electrons: 



— Wc(k&,P — a;k'a,$ — a) 

— a „ ^ 

kk'pp' 



(k 



x %+p-,fe'+p-'^( e ^ + e P-<y - e k'a ~ £p'- a )- (38) 



as well as equation of motion Eq. (|30p and resistivity 
tensor Eq. (|43|) derived before through phenomenological 
approach. 

Let us now describe briefly the idea of spin-grating 
experiments^ 9 -, where the spin-Coulomb drag has been 
observed. A periodic spin density can be induced by 
letting two linearly polarized light beams coming from 
different directions interfere on the surface of a two- 
dimensional electron gas. This interference produces a 
spatially varying pattern of polarization, with alternat- 
ing regions of left-handed and right-handed circular po- 
larization separated by linearly polarized regions. The 
spin density is optically induced in the regions of circular 
polarization. More precisely, the regions of right-handed 
circular polarization have a larger up-spin density while 
the regions of left-handed circular polarization have a 
larger down-spin density. At a given time t = the pump 
light is turned off and the subsequent time evolution of 
the spin-density is monitored by Kerr spectroscopy. In 
Kerr spectroscopy, one measures the amplitude of the 
spin-density modulation by looking at the rotation of 
the plane of polarization of the light diffracted by the 
spin-grating. The initial rate of decay of the spin grating 
amplitude 7 9 depends on the wave vector q of the grating 
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in the following manner^ 



7<z 



1 



D s q 2 



(39) 



where r s is the spin density relaxation time and D s is 
the spin diffusion constant. Therefore D s can be found 
from the slope of 7 9 vs. q 2 . The spin diffusion constant 
in the presence of spin-Coulomb drag was discussed in 
detail i n 64 ' 67 and has the following form: 



X°s/Xs 
1 + 7T 



(40) 



where D° s is the spin diffusion constant for the 
non-interacting system, x° s IS the spin-susceptibility 
for non-interacting system and \s is the interacting 
spin-susceptibility, with Landau-Fcrmi-liquid corrections 
taken into account. 

Actually, from the analysis presented in Refsi 64 i 67 one 
expects that the experimentally determined D s should 
include two effects: the Fermi liquid correction to the 
spin susceptibility 81 and the spin-Coulomb drag correc- 
tion. However, the Fermi liquid correction to the spin 
susceptibility is quite small. It is given by the well-known 
formula 8 -^ 



v° 

AS 

Xs 



1 



(41) 



where m/m e is the many-body mass enhancement and 
Fq is the Landau parameter described in detail in Ref4I. 
Since the m/m e ~ 0.96 and Fq ~ —0.2, the inter- 
acting spin susceptibility will be enhanced by no more 
than 20-30% and obviously will be independent of mobil- 
ity of 2DEG. Therefore, the Fermi liquid corrections to 
the spin-conductivity are very small in comparison with 
the spin-Coulomb drag corrections and the spin-Coulomb 
drag will be the main effect influencing the spin trans- 
port. Indeed, the experimentally determined Dj£ was 
found to be in excellent agreement with the theoretically 
predicted values for a strictly two-dimensional electron 
gas in the random phase approximatio n 64 ! 66 . Following 
this, Badalyan et aL 68 noticed that the inclusion of the 
finite thickness of the two-dimensional electron gas in the 



GaAs quantum well would worsen the agreement between 
theory and experiment, because the form factor associ- 
ated with the finite thickness of the quantum well reduces 
the effective electron-electron interaction at momentum 
transfers of the order of the Fermi momentum, which are 
the most relevant for spin Coulomb drag. Fortunately, 
it turned out that this reduction is compensated by the 
inclusion of many-body effects beyond the random phase 
approximation, namely local-field effects which, to a cer- 
tain extent, strengthen the effective Coulomb interaction 
by reducing the electrostatic screening. 68 The final up- 
shot of the more careful analysis is that the theory re- 
mains in quantitative agreement with experiment in a 
broad range of temperatures. 
IV. INFLUENCE OF SPIN-COULOMB DRAG 
ON THE EXTRINSIC SPIN-HALL EFFECT. 

A. Resistance tensor 

In this Section we study the influence of electron- 
electron interactions on the spin-Hall effect. Main discus- 
sion concerns 2DEG, however at the end of this Section 
we will comment on the behavior of spin-Hall conduc- 
tivity in bulk materials. We start from the Hamiltonian 
which includes electron-electron interactions: 



H — H n 



1 



2 4r*. e b \fi 



a 



£>x Vi^-* (42) 



where H ni is defined by Eq. and V ee = J2i^j e^-fjl ■ 
Notice that the electric potential coming from electron- 
electron interactions, like every potential whose gradient 
is non-zero, generates the spin-orbit term in Hamiltonian. 
This new spin-orbit term, introduces a new contribution 
to the Coulomb collision integral: 2eaha(E cr + E- a ) x 
z(k — k') which adds up to previous term i.e. the differ- 
ence of velocities for spin-up and down (see Eq. (j37|) ). As 
a consequence, electron traveling say in x direction with 
spin up can be scattered in a y direction with simultane- 
ous spin-flip, i.e. resistivity tensor contains terms which 
connect y and x components with opposite spins. The 
full resistivity matrix in the basis of x-\, y-f, xj_, y[ has 
the following form: 



/ pf + p bu njn ] 



Pf + A T pf 
P SD ~P'l 

7a 



D 



P\ - A TP| 

A \ P\ + P n ll n 1 

-Bj a 
-P SD P\ 



AT ~p SD 



p\ 

p D j_ „SD 



Br 
- p sd 



p-n r /n L -pf + \ lP f-A} 



•yet 



(43) 



Ala n D 
A l Pi 



+ p hU 7i ] /n i J 



where p SD — mf/ne 2 is the spin Coulomb drag re- 
sistivity and p' a = m/n^e 2 ^ (recall that X a = 
is a dimensionless quantity). AT and BJ a represent 



the terms of the first order in electron-electron cou- 
pling 7 and in SO coupling a and are defined as follows: 
AT = -\aPsDn-a/n a + 2maj[-n- a p® /n + (n_ ff /n - 
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n 2 _ a lnn a )p SD \ and BJ a = \ a p S D + 2ma-f[-n_ IJ p® cr /n + 
(n-cr/n — n a /n)p SD ]. Notice that the resistivity satisfies 
the following symmetry relations: 



13/3' 
Pa-a 



r — aa 



(44) 



(45) 



where upper indices (3 and /?' denote directions, and the 
lower ones spin orientations. New features of the resistiv- 
ity matrix|43]are the 7a and j 2 a terms, which appear in 
the transverse elements of the resistivity when the system 
is spin-polarized. Furthermore, the off-diagonal resistiv- 
ity elements p%y a are generally non zero. In the para- 
magnetic case (zero spin polarization) the aj 2 terms are 
zero and the resistivity matrix simplifies significantly. In 
this case, we find simple interrelations between currents 
and electric fields in the spin and charge channels. Omit- 
ting spin-flip processes (1/r' = 0) we obtain 



E c = P D 1 



E s = 4(p SD + p u )j s + 2(p ss - \p u - Xp SD )j. 



2(p ss - Xp D - Xp SD )j s 



X z 



D 



x z 



(46) 



(47) 



where the charge/spin components of the electric field 

are defined as E c — L ) E s = e\ — E^, and the 

charge and spin currents are j c = fa + j± and j s — JT ~ Ji , 
respectively. The spin-Coulomb drag renormalizes the 
longitudinal resistivity only in the spin channel. This is 
a consequence of the fact that the net force exerted by 
spin-up electrons on spin-down electrons is proportional 
to the difference of their drift velocities, i.e. to the spin 
current. Additionally, the electron-electron corrections 
to the spin-orbit interactions renormalize the transverse 
resistivity in the charge and spin channels, so the On- 
sager relations between spin and charge channels hold. 
Under the assumption that the electric field is in the x 
direction and has the same value for spin up and spin 
down electrons we see that Eq. (|4"6"|) and Eq. (|4"Tj) yield 
the following formula for the spin current jf y = j-\ — jy 
in y direction: 



■>s,y 



p ss /(p d ) 2 A 

l + p SD/ p D pD 



(48) 



The first term in the square brackets is associated with 
the skew-scattering, while the second is the side-jump 
contribution. Notice that the side-jump conductivity 
o B i = — = — 2ane 2 depends neither on the strength of 
disorder nor on the strength of the electron-electron inter- 
action. Moreover, as we showed ir*21 by using a gauge in- 
variance condition, the side jump does not depend on the 
electron-impurity and electron-electron scattering poten- 
tial to all orders in both these interactions. By contrast, 
the skew scattering contribution to the spin conductiv- 
ity, in the absence of e-e interactions scales with trans- 
port scattering time. Therefore for very clean samples, 
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FIG. 8: Spin Hall conductivity as a function of temperature 
at constant mobility. o-^j.(open black squares), ay% (open 
blue and close red circles) ,<t*°' (red solid line and dashed 
blue line) are the side-jump, skew-scattering and the total 
spin conductivity in a presence of electron-electron interac- 
tions. Notice that only skew scattering conductivity is mod- 
ified by spin-Coulomb drag. Red and blue curves/symbols 
are for /1 = lm 2 /V.s, and p = 0.1m 2 /V.s. We choose the 
parameters typical for the experimental 2DEG confined in 
Alo.1Gao.9As quantum well i.e. density of electrons and im- 
purities n2D = rii — 2.0 x 10 12 cm -2 , m — 0.074m e , and 
two sets of mobilities and relaxation times: p = 0.1m 2 /V.s, 
r = 4 x 10^ 5 ns t ss = 0.02ns and p = lm 2 /V.s r = 4 x 10 -4 
T sa = 0.2ns. The effective spin-orbit coupling ah = 4.4A 
in accordance with^i. We used the model potential (see 
appendix) where an effective impurity radius a — 9.45nm, 
the height of attractive impurity potential Vb = — 5meV for 
p = 0.1m 2 /V.s and F = -1.6meV for p = lm 2 /V.s. 



the skew scattering contribution would tend to infinity. 
However, this unphysical behavior is cured by the pres- 
ence of spin-Coulomb drag, which sets an upper limit to 
the spin conductivity of the electron gas: so the skew 
scattering term scales as r/(l + psd/pd) = t/(1 +7t), 
which tends to a finite limit for r — > 00. Let us now make 
an estimate of skew scattering contribution to the spin- 
Hall condutivity. Using the typical ratio of t/t ss w 10~ 3 
(see Eq. (|24p we obtains: 



10 



-3 a D 



1 + 7T 



(49) 



We will use this estimate in Section VI to compare im- 
portance of different contributions to the spin-Hall con- 
ductivity. 

Also, the direct dependence of the skew scattering con- 
ductivity (see Eq. (|49)l ) on the transport scattering time r 
is the reason why this term is modified by spin-Coulomb 
drag. By contrast, the side-jump conductivity, which is 
independent of r, remains completely unaffected. The 
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total spin Hall conductivity may either decrease or in- 
crease as a result of the spin Coulomb drag, depending on 
the relative sign and size of the skew-scattering and side- 
jump contributions. In the ordinary case of attractive 
impurities, when the skew-scattering contribution dom- 
inates, we expect an overall reduction in the absolute 
value of the spin Hall conductivity. 

This is shown in Fig. (8). One can see that in high- 
mobility samples the spin Coulomb drag reduces the spin- 
Hall conductivity very effectively. There is no upper limit 
to the reduction of the spin-Hall conductivity since the 
factor I + 7T can become arbitrarily large with increasing 
mobility. The behavior of jt as a function of temperature 
is shown in Fig. (9) for a typical semiconductor mobility, 
fi = 3 x 10 4 cm 2 /V.s. For example, in a two-dimensional 
GaAs quantum well at a density n — 10 11 cm and 
mobility 3 x 10 4 cm 2 /V.s the factor jt is quite signifi- 
cantly larger than one in a wide range of temperatures 
from 50 K up to room temperature and above, and can 
substantially reduce the skew scattering term. 
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FIG. 9: The ratio of spin-Coulomb drag to Drude resistance 
7T as a function of temperature for GaAs and InAs 2DEG 
with mobility fi — 3 x 10 4 cm 2 /V.s. The dashed and solid lines 



correspond to n = 10 



and n — 10 



respectively. 



Let us finally comment on the spin-Hall conductivity 
in 3D. Although the general formulas are the same as in 
2D, the actual value of p ss must be obtained by solving 
a three-dimensional scattering problem. This has been 
done in Refi 7 ^ for the following model attractive potential 
between the electron and an impurity atom: 



V(r) = 



(50) 



where e is permittivity of material and 1 / q s is the screen- 
ing length associated with the Thomas-Fermi screening. 



For this model potential, the spin-Hall conductivity takes 
the form: 



a SH = -2ane 2 + ^-a D 



(51) 



where ctd is the Drude conductivity and 7 S is a skewdness 
parameter, given approximately by 7 S = 4a/*i/(a 2 3 ). 70 
Notice that the side-jump contribution has exactly the 
same form as in 2D. Further, in 3D the spin-Coulomb 
drag should modify the spin-Hall conductivity in a simi- 
lar way i.e. by renormalizing the skew-scattering by fac- 
tor I+7T while leaving the side jump contribution to the 
spin-conductivity unchanged. Therefore except for the 
different scaling of 7 with temperature the spin-Hall con- 
ductivity behaves very similarly in 2D and in 3D . 65 i 66 



B. Spin accumulation 

A quantitative theory of the spin accumulation in semi- 
conductors requires in general a proper treatment of the 
boundary conditions as well as electron-hole recombina- 
tion effects^ 2 -. In this chapter we will study the influence 
of spin-Coulomb drag on spin accumulation assuming 
that electrons are the only carriers involved in transport. 
Our goal is to interpret the optical experiments in which 
spin accumulation is measured ( 39 i 40 ) . Notice that in pre- 
vious theoretical paper o 70 i 71 , directions of electric field 
and spin accumulation where exchanged in relation to 
experimental ones which led to a difference in a sign be- 
tween experiment and theoretical predictions due to the 

-o SH . In 

ryx 



following relation between resistivities: p^y 
this review we finally clarify this point and show that the 
sign of experimental and theoretical spin-accumulations 
agree. 

We consider a very long conductor in the form of a bar 
of length L in the y direction and narrow width W in the 
x direction, exactly the same setup as in experiments^ 4 ^ 
(see Fig (10b)). A charge current flows only in the y di- 
rection. The spin components of the transverse current 
j£, with a —\ or J, add up to zero everywhere and indi- 
vidually vanish on the edges of the system, i.e. j% = 
at x = ±W/2. In order to satisfy the boundary con- 
ditions the system cannot remain homogeneous in the 
x-direction. A position-dependent spin density, known 
as spin accumulation develops across the bar, and is re- 
flected in non-uniform chemical potentials fi a (x). In the 
steady state regime the spatial derivative of the spin- 
current in the x-direction must exactly balance the re- 
laxation of the spin density due to spin- flip processes i.e.: 



e dJ s _ p a {x) - p-g(x) 
a, dx LI 



(52) 



where L s is the spin diffusion length and a s is the longi- 
tudinal spin-conductivity. Additionally, Ohm's law must 
be fulfilled: 



J s = a s (E* - E*_ a ) 



(53) 
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where the effective electric field in the x- direction is 
equivalent to the gradient of chemical potential: 



(54) 



Notice that in the limit of infinite spin-relaxation time 
(L s — > oo ) the divergence of spin-current equals zero and 
the spin accumulation can be obtained directly from the 
homogeneous formulas, Eqs. (|46f and (I47)) . In an inho- 
mogenous case, combined Eq. (f52"j) and Eq. (f5*5|) lead to 
the following equation for the spin accumulations 83 



taken into account. This effect may increase or decrease 
the total spin accumulation, depending on the relative 
sign and magnitude of the side-jump and skew scatter- 
ing contributions. It reduces it in the common case, for 
an attractive electron-impurity potential, where the two 
contributions have opposite signs and the skew-scattering 
dominates. 



d 2 \}i a (x) - H- a (x)] p a (x) - Li_ a (x) 



d 2 x 

whose solution is: 



p<j(x) — P-a{x) — Ce L ° + C'e 



(55) 



(56) 



and C, C are constants to be determined by the bound- 
ary conditions j± cr (±W/2) = 0. Additionally using 
Eq. (j54"|) and the resistivity tensor we can write the 
boundary conditions for E%(±W/2). Using the bound- 
ary conditions for the spin-dependent chemical potentials 
and the spin-dependent electric fields one finally finds the 
following formula for the spin accumulation in a param- 
agnetic case: 



2eL s E y [p ss - Xp D - Xpsp] sinh(x/£ s 
p D cosh{W/2L s ) 



(57) 



The formula for the spin-accumulation in a spin-polarized 
case can be also easily obtained and the interested reader 
can find it in RefiZI. Finally, the spin-accumulation at 
the edges of sample for L = W/2 has the form: 



V: c = p^W/2L s )~p l (W/2L s ) = 
2eL sJy [p ss - Xp D - Xpsn] tanh(W/2L a ) 



(58) 



The three terms in the square brackets of Eq. (I58|) are 
the skew-scattering term, the ordinary side-jump contri- 
bution, and a Coulomb correction which has its origin in 
the side-jump effect. The latter is not a spin Coulomb 
drag correction in the proper sense, for in this case the 
transverse spin current, and hence the relative drift ve- 
locity of the electrons, is zero. What happens here is that 
the spin Hall current is canceled by an oppositely directed 
spin current, which is driven by the gradient of the spin 
chemical potential. Now the spin Hall current contains a 
universal contribution, the side-jump term, which is not 
affected by Coulomb interaction, but at the same time 
the constant of proportionality between the spin current 
and the gradient of the spin chemical potential, that is to 
say the longitudinal spin conductivity, is reduced by the 
Coulomb interaction. Therefore, in order to maintain 
the balance against the unchanging side-jump current, 
the absolute value of the gradient of the spin chemical 
potential must increase when the Coulomb interaction is 
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FIG. 10: (color online) (a) Spin accumulation as observed 
in experiments by Kato et al.— (b) the spin accumulation 
predicted theoretically (see Eq. ([58])). For experimental part 
of Figure: part (C) is a Kerr rotation as a function of x and 
external magnetic field for electric field E — WmVm , parts 
(D and E) describe spatial dependence of Kerr rotation peak 
Ao and spin lifetime r s across the channel, respectively. Part 
(F) shows the reflectivity R as a function of x. 



Additionally, Coulomb interactions affect the spin ac- 
cumulation indirectly through the spin diffusion length 
as shown in the equation below: 



L„ = 



Xs 1 + PSD I PL 



(59) 



which follows immediately from Eq. (j40|) . However, in the 
limit of W L s , tanh(TF r /2i s ) can be approximated by 
W/2L S , and the spin accumulation at the edges becomes 
independent of L s . In this limit, the influence of the 
Coulomb interaction on the spin accumulation is only 
through the XpsD term. 

Let us now put in some numbers. For a two- 
dimensional electron gas in an Alo.1Gao.9As quantum 
well^ with electron and impurity concentrations rii = 
n 2 D = 2 x 10 12 cm~ 2 , mobility ^i=0.1m 2 /V.s, L s = 1/xm, 
t = 4 x l(T 5 ns, t ss = 0.02ns, ah = AAA 2 , j x = 0.02 
A/cm and for the sample with width W = 100//m, we 
calculate the spin accumulation to be — 1.5meV/|e| on 
the right edge of the sample (relative to the direction 
of the electric field) i.e. for W/2 = 50/im. This means 
that the non-equilibrium spin-density points down on the 
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right edge of the sample and up on the left edge exactly 
like in the experiment. 

The inhomogenous profile of spin-accumulation is pre- 
sented in Fig. (10). Fig. (10a) shows the signal of the 
spin-accumulation (actually the Kerr rotation angle) ob- 
served in the experiment, while Fig. (10b) shows the 
profile of spin- accumulation expected from the formula 
1551 The general profile of the spin accumulation is 
satisfactory on a qualitative level (taking into account 
that we considered a very simple description of spin- 
accumulation) . 




0,01 -4-r™, 

0,1 1 10 100 1000 



T[K] 

FIG. 11: (color online) Mobility, fi, as a function of tem- 
perature, T, for three different low-T /i's. In inset, the spin 
accumulation (Va C ) vs. T. The side jump contribution to V^ c 
dominates for low T. For increasing T, the lower temperature 
red cross corresponds to the T where the sign of V^ c starts to 
be controlled by skew scattering, the higher temperature red 
cross to the place where side jump dominates again. 

As we mentioned before, it is possible to distinguish 
between side jump and skew scattering contributions 
to the spin accumulation because they scale differently 
with mobility. We have proposed an experiment to dis- 
tinguish between these two contributions in a study of 
the temperature dependence of o**3 or spin accumula- 
tion V* c = —Vaci where the last equality stems from 
the fact that the transverse resistivity elements are odd 
under exchange of spatial coordinates i.e. p xy = —p yx . 
Fig. (11) presents the behavior of mobility versus tem- 
perature for experimentally attainable samples. Due to 
different scattering mechanisms, the mobility scales non- 
monotonically with the temperature. Hence /x will grow 
as T 3 / 2 for low T as a result of scattering from ion- 
ized impurities and will decrease as T~ 3 / 2 for larger T 
due to phonon scattering. It is thus possible to observe 
two changes of sign of V£ c moving from low to high Ts. 
fj, = l/(AT~ 3 / 2 + BT 3 / 2 ), where A was found from the 
low-T mobility and B was fixed by a room temperature 
mobility of 0.3 m 2 /Vs for AlGaAs. At low T the mobility 
is low and the side jump contribution to V y c dominates. 



With increasing T, the first cross designates the point 
where the skew scattering begins to dominate, and the 
second cross, at higher T, is the point where the side 
jump takes control of the sign of V£ c again. Even if the 
sign change is not detected, by measuring whether V y c 
increases or decreases as \i increases with changing T it 
should possible to tell whether side jump or skew scatter- 
ing dominates. Notice that the values of parameters for 
the theoretical curve designated by circles are exactly the 
same as the values reported for the samples in the recent 
experiments^ - on a [110] QW. The samples with lower 
mobilities can be easily obtained by additional doping 
with Si inside the quantum well. 

V. INFLUENCE OF RASHBA TYPE 
SPIN-ORBIT INTERACTION ON 
SPIN-COULOMB DRAG 

So far, the influence of the Coulomb interaction on 
the intrinsic spin Hall effect has not been analyzed. The 
issue is more complex than the problem presented in Sec- 
tion IV i.e. the study of the influence of spin-Coulomb 
drag on the extrinsic spin Hall effect. The main differ- 
ence between the problem with impurities and the prob- 
lem that considers spin-orbit interaction coming from the 
band structure is that the latter usually does not conserve 
the z-component of spin. The general form of intrinsic 
spin-orbit interactions in 2D is: 

H b = -h(k)-a (60) 

where a is the vector of the Pauli matrices and b(k) is the 
intrinsic spin-orbit field. Due to time reversal symmetry 
this field needs to fulfill the following condition: b(k) = 
—b(—k). For example for the simplest model describing 
the spin-orbit interactions in a 2DEG oriented in [001] 
direction one has^ 8 -: 

b(k) = 2az x k (61) 

where z is the unit vector in z direction and a is the 
spin-orbit coupling strength. This model is known in lit- 
erature as the Rashba model^ 8 -. For 2DHG, the Luttinger 
Hamiltonian for low densities can be simplified by taking 
into account only the heavy-hole ban d 31 ! 84 . This gives 

H 2 DHG=ia h {k 3 _a + -k\a-). (62) 

Obviously, these Hamiltonians do not conserve the z- 
component of spin and therefore the simple approach pre- 
sented in Section IV can not be applied. The main com- 
plication is that we need to consider the whole density 
matrix (also the off-diagonal terms) in the Boltzmann 
approach. However, the Rashba model (at least for the 
standard definition of spin-current^) gives zero spin-Hall 
conductivit y 29 ' 30 ! 32 ' 33 ! 34 ' 38 so we do not expect that this 
result will be further modified by e-e interactions. On the 
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other hand, for the cubic spin-orbit models, vertex cor- 
rections are not important and the spin-Hall conductivity 
is of the order of e 2 jh. In this case, our expectation is 
that the spin-Hall conductivity, being non-universal, will 
be reduced by spin-Coulomb drag. However, there has 
not been yet the calculations which would quantitatively 
address this problem. 

So far the only studied problem was the influence of 
spin-orbit interactions on the spin-Coulomb drag 8 -^. Dia- 
grammatic calculations^, suggest that the spin-Coulomb 
drag is actually enhanced by spin-orbit interaction of 
the Rashba type, at least in the weak scattering regime. 
The correction is simply additive to what we would ex- 
pect from microscopic calculations of spin- Coulomb drag 
without spin-orbit interaction^ 2 -. It has the form 3(7* nt ) 2 
where j* nt = am e / {vftti) where vf is the Fermi velocity. 
There is still much room of course for more detailed stud- 
ies of the interplay between spin-orbit interaction and 
spin Coulomb drag. Recently, Weber et alM- have un- 
dertaken experimental studies of the relaxation of a spin 
grating in a two-dimensional electron gas oriented in such 
a way that the spin-orbit interaction is relevant. A theo- 
retical study by Weng et al.~- which takes into account 
both spin-orbit coupling and electron-electron interaction 
concludes that the spin Coulomb drag will still be visible 
as a reduction of the spin diffusion constant. The latter 
is still determined from the initial decay rate of the am- 
plitude of the spin grating, but the full time evolution of 
the amplitude involves two different spin relaxation times 
for in-plane and out-of-plane dynamics respectively. 



VI. EVOLUTION OF SPIN-HALL EFFECT 

As we showed in previous Sections, the spin-Hall ef- 
fect has various contributions. Therefore, it is in place 
to compare the importance of different mechanisms con- 
tributing to the spin-Hall conductivity as a function of 
hr/m. We chose hr/m as a SHE evolution parameter for 
two reasons (1) it has a dimension of a squared length and 
can be directly compared with the strength of the spin- 
orbit coupling ah (2) hr/m can be easily connected with 
the mobility. When hr/m is expressed in A 2 it is approx- 
imately equal 6/i, where /i is the mobility in cm 2 /V.s. 

In a d.c. limit we can distinguish three different 
regimes 89 : (1) ultraclean regime where - <C E so , where 
E so is the spin-orbit energy scale defined as E so — 
Ep{ah/ 'a 2 B ), as is the effective Bohr radius (2) clean 
regime characterized by inequality E so <C ^ -C Ep, 
where Ep is the Fermi energy, (3) the dirty regime in 
which - > Ep. In terms of hr/m these three regimes cor- 
respond to (1) £ » 2* (ultraclean), (2) a% < « £f 
(clean) and (3) — < a 2 B (dirty) where we assumed 
n = a^ 3 - On top of these limits we need to know the 
relative importance of the skew scattering and side jump 
contributions to the spin-Hall conductivity. Let us there- 



fore recall the final formulas for the skew scattering: 
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-3 °D 



1 + 7T 

and side jump spin-Hall conductivities: 



-2ane 



The ratio of these two conductivities is: 

a sj 2 ah\A 2 ] 

a ss \i\cmr I V.s\ 



(63) 



(64) 



(65) 



where we used the connection between mobility and 
hr/m. Also, in the above formula ah is in units of 
A 2 while fi is in cm 2 /V.s. The side jump and skew 
scattering contributions have the same magnitude when: 
— = 2 x 10 3 a?i. As we already mentioned the skew scat- 
tering contribution scales with the mobility and therefore 
will dominate the ultraclean regime, where this quan- 
tity is the largest. The skew scattering contribution is 
cut-off by spin-Coulomb drag when 7 3> 1/r and the 
cut-off value of spin-conductivity is cr^u* = 10 In- 
trinsic contribution (if it is not zero, like in a Rashba 
model) dominates in a clean regime. As was shown 
theoretically 9 -^, the vertex corrections connected with dis- 
order are zero for p-doped semiconductors and the spin- 
Hall conductivity is of the order of e 2 /h and therefore 
much larger than the side-jump contribution. For ex- 
ample for typical 2D semiconducting hole gases with 
densities (10 cm ), side jump contribution is around 
thousand times smaller than intrinsic one. Therefore 
in the scenario where the skew-scattering regime passes 
to the side jump regime the scale 2 x 10 3 a?L must be 

a 4 

larger than --§■. For doped semiconductors, for exam- 
ple GaAs, the effective Bohr radius is 100A, so we have 

|«4x 10 6 a?i > 2 x 10 3 a?i. Therefore we will have a 
direct transition from the skew scattering to the intrin- 
sic regime. However, there are two options to observe the 
side-jump effect: (i) the intrinsic contribution is zero (like 
in Rashba model), then side-jump could be observed in 
the clean limit (ii) the mobility is decreased by changing 
the temperature and the side-jump contribution could be 
observed in ultraclean regime as we mentioned in Section 
V (see Fig. (11)). In the dirty regime the spin-Hall con- 
ductivity diminishes to zero. The scenario of evolution of 
spin-Hall conductivity in semiconductors with non-zero 
intrinsic contribution is presented in Fig. (12). In con- 
trast, if one adopts Eq. (|6"3")) and Eq. (|64p to describe 
the spin-Hall effect in metals, one finds that the skew 
scattering term would evolve into side-jump contribution 
for hr/m* = 10 3 a and the intrinsic effect will eventu- 
ally appear for hr/m* = a B /(ah) in the clean regime. 
For parameters typical for Pt, the side jump and intrin- 
sic contributions are of the same order and dominant. 
Therefore further calculations (including the complexity 
of band structure) and experiments are needed to distin- 
guish between various mechanisms contributing to SHE 
in metals. 
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FIG. 12: (color online) Evolution of the spin-Hall effect as a 
function of Tit /m assuming that an electron-impurity scatter- 
ing potential has the form of the screened Coulomb potential. 
The side jump could be observed in an ultraclean regime if 
we decrease the mobility by for example the tuning of tem- 
perature. 



feet and the spin Coulomb drag effect, and the relation 
between them. Careful readers will notice that there are 
still plenty of open questions and unsolved problems. 



From the theoretical point of view, perhaps the most 
urgent open challenge is the calculation of the influence 
of the spin-Coulomb drag on the intrinsic spin Hall effect. 
From the experimental point of view, it would be inter- 
esting to see time-resolved studies of the spin Hall effect, 
possibly conducted by spin-grating techniques'^ 7 or by 
optical spin injection technique s 43 ' 44 . Furthermore, a di- 
rect detection of the influence of the spin- Coulomb drag 
on the spin-Hall effect and an experimental verification 
of the theoretical predictions of Section IV would be of 
great interest. Finally, a full description of the interplay 
between spin-orbit coupling and spin Coulomb drag re- 
mains an open challenge, particularly at the experimental 
level. 



VII. SUMMARY 



In this review we have summarized the current status 
of the knowledge concerning the extrinsic spin Hall ef- 
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